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ABSTRACT 


The amount of energy absorbed when a stream of external radiation enters 
a rectangular-groove cavity has been calculated for a variety of surface 
conditions, cavity depths, and incident energy distributions <> The surfaces 
of the cavity are either diffuse reflectors or specular reflectors, and a wide 
range of reflectivity values have been considered » The incoming radiation 
is either diffusely distributed across the cavity opening or else, arrives in 
a buhdle of parallel rays, As a function of these parameters, results are 
reported in terns of an apparent absorptivity, which is the ratio of the 
energy absorbed in the cavity to that which enters. The results show that 
for diffuse incoming radiation, a specular cavity absorbs more effectively 
than does a diffuse cavity. For incoming radiation in a parallel ray bundle, 
the comparison depends on the angle of inclination of the rays, 

NOMENCLATURE 

A surface area 

B rate of radiant flux leaving unit surface area, radiosity 
D, E constants of integration 

rate of diffuse radiation per unit area of opening 
F angle factor 

0, H constants of integration 

H rate of incident energy per unit surface area 
h cavity width 

numerical constants, eq, (18) 
k number of surface contacts 

L cavity depth 

directly “illuminated depth 
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exponents 

overall rate of energy absorbed 

local rate of energy absorbed per unit area 

energy rate per unit area normal to ray 

equivalent temperature of black surface, ~ 

dimensionless coordinate, x/L 

distance from cavity opening 

separation point for k and k~ 1 surface contacts 

dimensionless coordinate, y/L 

distance from cavity opening 

surface absorptivity 

apparent absorptivity of cavity 

dimensionless radiosity, B/e^ or B/S sin y 

inclination angle, fig. 1 

numerical constants, eq„ (18) 
distance, Jx ■» Y| 
reflectivity, 1 ~ u. 

Stefan -Boltzmann constant 



INTRODUCTION 


This paper is concerned with the energy absorbed when radiation from an 
external source enters a rectangular-groove cavity. It is well-known that the 
energy absorbed within any cavity is greater than that which would be absorbed 
by a plane area stretched tightly over the opening of the cavity* This is called 
cavity effect and may be explained in terms of the additional opportunities for 
energy absorption which accompany the multi-reflections within the cavity. The 
quantitative extent of the cavity effect depends on the shape of the cavity,, 
as well as on the reflecting properties of the surface and on the nature of the 
incoming energy. 

The aim of this investigation is to quantitatively determine the cavity 
effect for the rectangular groove. The analysis is extended over a wide range 
of conditions which include surface reflectivity properties, nature of incoming 
radiation, and groove aspect ratio. Consideration is given to surfaces which 
may reflect either diffusely (i.e., according to Lambert’s Cosine Law) or 
specularly (angle of incidence equals angle of reflection). For each one of 
these reflectivity characteristics, consideration is given to radiation which 
enters the cavity either as a diffuse stream or as a bundle of parallel rays. 

For these various cases, energy absorption results are obtained and repotted as 
a function of surface absorptivity, groove aspect ratio, and inclination ; ngle 
of the parallel ray bundle. 

A schematic diagram of the rectangular-groove cavity is shown in figure 1 
along with dimensional nomenclature. The depth of the groove is measured by L, 
while the width is h. The length of the cavity in the direction normal to the 
plane of the figure is sufficiently great so that end effects are negligible. 

In the case of a parallel ray bundle entering the cavity, the angle f is used 
to specify the inclination of the rays to the normal. 

It will be convenient to give separate consideration to the cases where 
the incoming energy is diffusely distributed or arrives in a bundle of parallel 
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rayso The former will be treated first and the latter second* Separate 
treatment will also be given for diffusely-reflecting and specularly-reflecting 
surfaces* 


DIFFUSE INCOMING ENERGY, DIFFUSE SURFACES 


The stream of diffuse incoming energy is uniformly distributed over the 
opening of the cavity and may be characterized by e^ per unit area of the 
opening. For purposes of analysis, the condition of diffuse incoming energy is 
exactly enui valent to that provided by a plane black surface of emissive power 
e^ and temperature (e^ » rgT^) stretched tightly over the cavity opening 
and radiating into the cavity. 

Consideration is first given to the case where the walls of the cavity are 

diffusely reflecting* The first step in the analysis is to write a radiant flux 

balance at an area dA (see figure l) located at a typical position x * Setting 

x Q o 

aside the contribution of emission,* the energy leaving dA is simply the 

x o 

reflected portion of the incident energy. Denoting the energy leaving a surface 
location per unit time and area by the radiosity B, and the incident energy per 
unit time and area by H, we can write 

B(x ) * (1 - oc ) H(x ) (1) 

O O 


where «*. is the absorptivity and p - (1 -rx) is the reflectivity* The energy 

incident »t x arrives via two paths: (a) directly from the external source 

o 

and (b) indirectly due to energy r eflected at the other surfaces of the cavity* 


From the external source, there is incident per unit area at x^ 

e. F . 
i x 0 =i 


(2a) 


where F . is an angle factor representing the fraction of the energy leaving 
o 1 

Energy from an external source enters the cavity, and it is desired to know the 
amount which is ultimately absorbed* With this in mind, it is only necessary 
to consider the radiant interchange process for energy which originated at 
the external source, and energy emitted by the surface need not be considered* 
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dA which arrives at the surface i stretched across the cavity openings From 


a typical position x (arva dA^) on the surface of the cavity, there is incident 


per unit area at x 


B(x) dF 


x o“ x 


(2b) 


where dF is tie fraction of the energy leaving dA which arrives at dA 0 

X ■" X X 

o o 

However, dA receives radiation from all surface locations within its range of 
x Q 

visibility, a*d the total contribution is obtained by integrating expression (2b) „ 
Thus, 'ith equations (2a) and (2b) for the incident energy, the radiant flux 
balance (") becomes 


B(x o ) = (l~o0 


_ e i 


\B(x) 

■x 


dF. 


v x 


(3) 


By in ^ pact ion of this equation, it may be noted that the unknown B appears under 
the integral si ;n as well as in other parts of the equation.. Equation (3) is 
therefore an integral equation.. 

Once solutions for the radiosity B are available, the energy absorbed in the 
cavity can ci calculated., First, ^he rate q at which energy is locally absorbed 
(per unit .rea) is simply found by multiplying the incident energy H by the 


absorptive by, 


otH 


c*. 

1 ~o< 


B 


(U) 


vihire equation (l) has been used* Next, the total rate Q at which energy is 
.-.bsorbec in the cavity as a wnoie may be obtained by integrating equation (U) over 
the surface of the cavity <> For a unit length in the direction normal to figure 1, 

Tbdx <5a) 


where, due to symmetry, the integral is extended over only half the surface and 

a factor of two has been included., It is convenient to report the results in 

terns of an apparent absorptivity defined as 

a 

o< ~ (total absorbed energy)/ (total incoming energy) (6) 

St 

Noting that the energy entering the enclosure per unit length normal to the 
plane of figure 1 is eh, the apparent absorptivity may be evaluated utilizing 



aquation (5) as 


0< a = -T~5T j h p 




where » B/e^ and X = x/L» It is thus seen that the solution of the integral 

equation (3) holds the key to the determination of the apparent absorptivity® 

In approaching the solution of this integral equation, it is advantageous 

first to simplify its form® Noting that energy leaving dA must either strike 

x o 

the other surfaces of the cavity or else pass out through the opening, it follows 
that 


dF 


= 1 


Introducing this into the integral equation (3) and rearranging.* there is obtained 


B'(x ) * <xe. + (l-oO \ B* (x) dF. 

° 1 Jx 


X 0 ~X 


where 


- B, or p « 1 - (B* /e i ) 


(7) 

(7a) 


But, the integral equation for B* is precisely that which describes the emission 
problem for a cavity having gray, diffuse walls which are at a uniform temperature 
(e^ as o»T^)o Moreover, solutions for the emission problem ars already avail- 
able (reference 1), and these may be carried over to the absorption problem being 
considered here by application of the simple relationship indicated in equation (7a)® 
The emission solutions covered the range of cavity aspect ratios L/h from 
0®25 to 10 for surface absorptivity values o ( of 0®5, 0®75* and Results 

for the apparent absorptivity o( a corresponding to these conditions hare been 

calculated by application of equation (5a)® Utilizing these, o( values have been 

a 

graphically interpolated to good accuracy for c( = Oo7 ; , in order to provide com- 
parisons with later results of this paper® The apnarmt abscrptivities thus 
obtained have been plotted on figure 2 as a function f cavity aspect ratio for 
parametric values of the surface absorptivity® 

Inspection of the figure reveals that for cavities which are not too deep, 
the apparent absorptivity increases with increasing dipth of cavity® However, for 
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sufficiently deep cavities, the apparent absorptivity becomes independent of any 

further increases in depth. The depth of cavity required to achieve the linn. ting 

value of apparent absorptivity decreases as the surface absorptivity increases » 

For instance, f or o( = 0.9, the limiting o( value (0.976) is essentially ac.iiered 

a 

at L/h — 2.5} while f or o( - 0.5, a hole depth L/h = 7 is required to essentially 
achieve the limiting value of (0.850). 

These findings are made plausible by considering the nature of the diffuse 
reflection process. The energy density of diffusely reflected radiation is the 
same in all directions. Therefore, some of the energy reflected at each surface 
contact passes outward through the opening of the cavity, regardless of the depth 
of the cavity. It is for this reason that all of the energy entering the cavity 
cannot be fully absorbed. 

In appraising the results, it may be noted that in tht absence of the cavity 

effect, cx would be identically equal to ert . Thus, the increase of c< relative 
a a 

to oc is a measure of the magnitude of the cavity effect. Tea results cf figure 2 
indicate that the cavity effect is most pronounced for deep cavities and for 
surfaces of low absorptivity. Thus, utilizing a diffuse cavity having an aspect 
ratio L/h s- 7, a ?0£ increase in the effective absorbing power can be achieved 
relative to a plane surface of absorptivity 0.5<> Moreover, a 60?' increase oan 
be achieved with a cavity whose aspect ratio is only 1.82 

DIFFUSE INCOMING ENERGY, SPECULAR SURFACES 
Consideration is now extended to a cavity whose surfaces are specular 
reflectors on which the angle of reflection equals the angle of incidence. 

It will be convenient to replace the incoming diffuse radiation by a black surface 
of emissive power e^. With this, there is formed a co inlete enclosure composed 
of the three specular walls of the cavity (designated zs 1, 2, and 3) and the 
black surface (designated as U)« A schematic diagram of the enclosure is shown in 
figure 3c 



6 


One approach toward solving for the radiant interchange within the enclosure 
would be to keep account of individual rays as they reflect and re-reflect on the 
specular surfaces. Fortunately, this arduous task may be circumvented by applying 
a new formulation (reference 2} which utilizes a basic property of images which 
are formed in plane mirrors. Namely, that light or radiant energy reflected from 
a plane mirror appears to come from an image which is located behind the mirror, 
the distance between the image and the mirror being the same as the distance 
between the object and the mirror. 

The development of tne new calculation method is described in the reference, 
and only its application to the absorption problem in the rectangular cavity will 
be discussed here. Considering figure 3, it is easy to see that the energy absorbed 
at the cavity walls is equal to the difference between the energy emitted by the 
black surface and that which is absorbed by it. The energy emitted per unit 
length normal to the plane of the figure is e^h. Then, denoting by the energy 
incident per unit time and area on the black surface, it follows that the rate Q 
at which energy is absorbed in the cavity as a whole is 

Q = h 1 (8a) 

In terms of the apparent abeoip th r b ty o( defined by equation (6), this becomes 

5 . 

<*=l~(H./e) (8b) 

a lj i 

It is thus seen that the apparent absorptivity will be determined as soon at the 

incident energy H. has been calculated, 

u 

In the process of specular reflection by the mirrors 1, 2, and 3* a series of 
images of surface k are formed as indicated in figure 3° The notation may b e 
explained as follows: b?" is the image of surface U formed by one specular 

reflection, h is the image formed by two specular reflections, and in general, 
h n is the image formed after n specular reflections. The subscript u refers tc 
upper set of images, while the lower images are without subscripts. Clearly, 
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since the upper images lie in the same plane as surface 1*, energy leaving these 


images will not impinge directly on surface The radiation leaving a typical 


n 


lower image surface It per unit area is 

\ (1 - o() n 

where (l-oc) = p is the reflectivity.. Of this, an amount 

e (l-oc) n F (9a) 

U- it 

reaches surface it, where it is absorbed. The angle factor ^ corresponds to 
diffuse radiation between surfaces It and it n « From all of the image surfaces, 
there arrives at It — «*> 

e > (>-<*) F^,, t £[/' F ,.^ (9b) 

and ^ith this, the apparent absorptivity may be evaluated from equation (8b) as 

CO _ 

°<c>~ / — -i- 2>(/ h+.y (1 °^ 

L o-i ” - 

It only remains to furnish the angle factors F. , n . Wiese may be nmst easily 

it- it 

derived by the method of crossed and uncrossed strings due to Hottel* (reference 3). 
In this method, strings are imagined to be tightly stretched between the end 
points of surfaces it and lt n . The angle factor is then obtained by summing the 
lengths of the two crossed strings and subtracting away the lengths of the two 
uncrossed strings, and then dividing the difference by 2ho In this way, there 
is obtained 


it ~ If 


n 


i. 

2 


^ (n- 2) 2 + (2L/h) 2 +^j n 2 + (2L/h) 2 - 2 ||(n~ l) 2 + (2J./h) 2 

( 11 ) 


Inspection of equations (10) and (11 ) reveals that the apparent absorptivity 
OC depends upon two parameters, the cavity aspect ratio L/h and the surface 

absoiptivity Oi 0 For each h/L, cx. combination, the infinite series appearing 

in equation (11) was summed, and <X thus obtained. The actual numerical 

a 

* 

This method applies to surfaces which are very long in one direction. 
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calculations were performed on a Uni vac 1103 digital compute To The series was 

«,c 

truncated when trie contribution of a given term was less than >0 of the sum of 
the series. 

Results fcr the apparent absorptivity have been plotted on figure U as a 
function of cavity aspect ratio for parametric values of the surface absorptivity. 

In recognition of the fact that specular surfaces may have higher reflectivity 

✓ 

values than diffuse surfaces, lower c* values have been included on figure h° 

From the figure, it is seen that the apparent absorptivity increases monotonically 

with increasing cavity depth, ultimately approaching the value of unity for 

sufficiently deep cavities. This is in contrast to the behavior of diffuse 

surfaces (figure 2), for which the apparent absorptivity approaches a limiting 

valufc less than unity for very deep cavities. These contrasting findings may be 

understood by taking cognizance of the different nature of the two reflection 

processes. In the case of diffuse reflection, a portion of the energy reflected 

at nach surface contact passes outward through the cavity opening. In the case 

of specular reflection, only energy which reaches the bottom of the cavity and is 

ultimately 

reflected from the base surface 3 (figure 3) can ever/pass outward through the 
cavity opening. Thus, for sufficiently deep specular cavities, essentially all 
of the incoming radiation is absorbed. It is also easily understood that deeper 
cavities are needed to achieve the total absorption condition for surfaces of 
lower absorptivity. 

It is interesting to compare the absorbing power of diffuse and specular 
cavities for diffuse incoming radiation. Comparing corresponding curves on 
figures 2 ard it is seen that a specular cavity absorbs more effectively than 
does a diffuse cavity. The advantages of the specular cavity are accentuated for 
deep cavities and for low surface absorptivlties. For example, for o( o.5> the 
apparent absorptivity values for specular and diffuse surfaces at 1,/h - 1 ars» 
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0o? 62 and 0®7U2® At L/h = 10, the corresponding values are 0.963 and 0,850® 

A similar comparison shows smaller deviations for oi = 0, 9, 

PARALLEL INCOMING RAYS, DIFFUSE SURFACES 

When a parallel ray bundle enters the cavity under an inclination angle ^ , 
a length L^ of one cavity wall is directly illuminated as shewn in the inset of 
figure S» From the geometry of the figure, 

L^h/tan f (12) 

The remainder of the cavity walls receives radiation only be reflection® Different 
quantities of energy will thus be incident on oppositely-located positions on the 
parallel walls of the cavity® In cognizance of this lack of symmetry, separate 
coordinates x and y are employed to specify positions on the parallel walls, see 
figure So 

When the surfaces of the cavity are diffuse reflectors, the radiant inter- 
change withih the cavity is described by integral equations® Because of the 
a3symetry noted above, several integral equations are required® Additionally, 
tne solutions of these integral equations will depend on three independent para- 
meters: the aspect ratio L/h, the surface absorptivity c< , and the inclination 

angle y* of the incoming rays. A complete set of solutions to this problem involving 
several simultaneous integral ^equations and three independent parameters is a very 
formidable computational undertaking, especially since numerical means must be 
utilized. Now, in the previous case of diffuse incoming energy and diffuse sur- 
faces, it was found that cavities of moderate depth displayed absorption character- 
istics very little different from those of a cavity of infinite depth® Taking 
cognizance of this, it was decided to reduce the number of parameters in the 
present problem of incoming parallel ray3 by focusing attention on the cavity of 
infinite depth® Such a cavity is shown in the inset of figure 5® Approximate 
auxiliary calculations were carried out to determine the depth of the finite 
cavity for which the infinite cavity results would apply. 
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The mathematical formulation of the problem follows parallel to that which 
has beer, carried out for the case of diffuse incoming energy The radiant flux 
balance of equation (1) continues to apply. It remains to evaluate the incident 
energy H. If the energy carried by the parallel ray bundle is denoted by S per 
unit time and area normal to the ray, then 

S sin (13) 

is the rate at which external energy is incident per unit area of surface in the 
range 0 On all other regions of the cavity surface, the direct incidence 

of external radiation is zero. The contribution to H of radiation reflected 
from the surface of the cavity is obtained by integrating expressions having the 
form of equation (2b). Thus, we can write 



O < y^L, 


(Ilia) 


, (liib) 

L, oo 


Q <. <-jS. OO 


(Ihc) 


To complete the formulation, it is necessary to derive the angle factors a 

For the case where the surfaces are very long in the direction normal to tht plane 

of the figure, there is a particularly simple method for finding angle factors 

which is based on equation (31-58) of reference h° For instance, to find dF , 

x- y 

the normal is constructed at x and the connecting line is drawn from x to y. 

The angle included between the normal and the connecting line is denoted by <p . 
Then, applying the differential form of equation (31=58), which states that 

dF *- y = 1“ ^ d y 


there is obtained 


li 


dF 


y 


1 

2 r 


L( 


y- xf * h 2 


fj 3/2 


dy 


(15a) 


J 


A similar expression is found for dF , except that dy is replaced by dx. 

y ™ x - 

terms of dimensionless variables X - x/h, Y = y/h, the angle factors become 


In 


dF 


x » y 


[(T-X)2 + l]‘ 


dF 


y- x 


1/2 

|_(Y= X) 2 + l] ^ 


dX 


(15b) 


The function 1(1- X) 2 + IJ has been plotted as a solid curve on figure 6 S 
with ^ representing the distance Y-X . It is seen that this function has the 
general form of an exponential = The figure shows that a given surface location 
is most affected by radiation from other surface locations which are directly 
across from it* i 0 e 0 , at small ^ . 

A final dimensionless form of the integral equations (1U) is obtained by 
letting 

P « B/S sin £ (16) 

The effect of this substitution is to replace B by A in all terms of equations (lit), 
and to replace S sin ^ by 1. 

The numerical solution of equations (lit) as a function of two independent 


nq r*o mo fo r»c *S o r*+A 1 1 

f w afuw U VJ.AJ. 


formidable coniputatlonax uituv <> 


The lengthy digital 


computer solutions can, however, be circumvented by applying a method first used 
by Buckley (reference 5) with very good accuracy (reference 6) in connection with 
the emission problem in circular holes. Later investigators have used Buckley’s 
method for various cases of radiant interchange in circular-cylindrical tubes, 
but to the knowledge of the authors, this approach has yet to be applied to the 
rectangular groove. Buckley recognized that the angle factor varied approximately 
as an exponential function of distance, and he proposed that the actual angle 
factor be approximated either by a single exponential or by a sum of exponentials. 
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Following through on this idea for the case of the rectangular groove, one and 
two term exponential approximations for the angle factor have been derived and 
are plotted on figure 6» The two term exponential provides a good representation 
for the angle factor and will be used here. 

In incorporating the exponential approximation into the analysis, cognizance 
must be taken of the fact that]T represents a distance between points and is 
therefore intrinsically positive,. This means that for Y c X, £ denotes X- Y, 
while for Y > X, ^ denotes X» This point may be illustrated by evaluating 
equation (Ilia) in terms of the exponential angle fact or c 


ft(X) 


-"MX-yil ;u 

ic / y 


(i-a) +(i -«) | pc y;^K’ i £ “ ^ J ‘+Kie 

+(/-<*)( ( 3 <y) K,e' A ' iV ~ x l*&e Al<V X> d}/ 

)w= v 


where 


K »Xo2lt65 s lo203 

1 1 

X 2 * -Oo2u65 7 ^ « 6,8157 


(17) 


as: 


Integral equations similar to (17) can be written for aquations (lbb) and (llo), 
Continuing with the analysis „ it is convenient to represent the unknowns 
by four functions (b . * A • A « A , as follows 


Pi* (V (V fh 

P (X) « p X (X) , 

P (X) = |4g(X) , 
p (I) - (3 3 (Y) , 


co 


My) 


Pi 


(?) 


OiX^L./h 

i 

L ± /h < X £ 

0 4K L i A 
Lj/h ^ Y^od 


(19) 


The next step is to reduce the governing integral equations to differential 

or Y 

equations. This is possible because the X/dependenee of exponential, functions 
such as e" ^»( x " Y ) ^ e “ y ) is UJtc hanged by differentiation = In recog- 

nition of this, equation (17) is differentiated four times with respect to X. 


Then* by properly combining equation (i?) with the second and fourth derivatives : 


it is found that all the terns containing integrals cancel out, leaving 
- f A* -i- frl ) p>‘‘ + XXi p, - (i -&}( X, ), '• Kz 'Up?*' 

-(/ - (l -<*) /Chi 

where the primes denote differentiation. In a similar way, differentiation of 
equation (lire) with respect to I, Q I ^ L./h, yields 


(20) 


a / 'i 1 l /- '' 

~( A. t A 3 .) 


At /% - (/ " j( k, -t Kz ) iX j p> ( 


/ » / l \ ^ \ vs 

'(/ ~ Oi! fU A, & . 


o 


It is seen that equations (20) and (21) are a pair of s imultaneous , linear, 
fourth -order, ordinary differential equations for the functions and ^ - The 

ft 

homogeneous part of such equations must have exponential solutions 

D eP T 




ft 


( 22 ) 


Ini, reducing these into the homogeneous part of the differential equations (20) 
and (21) yields two linear, homogeneous , algebraic equations for S and I with 
coefficients which contain p, ) , Kp and A non- trivial solution to 

these homogeneous algebraic equations is possible only if the determinant of the 
coefficients is zero. This yields a polynomial of eight degree (the- so-called 
secular equation) for the exponent p, for which it is found that there are eight 
real roots. These roots pair off together as follows? + p^, + p^, + p^ 9 + 

A listing of the expressions for p, , p^, Py and p^ is given in the Appendix. 
Returning to the linear algebraic equations, it is additionally found that 

D ± = E , for i a 1, 2, 3, U 
\ » -E i , for i a 5, 6, 7, 8 

Therefore, with this information, the solutions for ft and A. can be written. 

i 1 j a 


right-hand sides equal zero. 


Ik 


3 _ ~ E + E e 1 ' 1 * + S 0 e^ X + E 
1 X O 1 2 ; 


-RX 


+ E. e‘^ X + E_, + 5, 


+ E + Eg 


(22a) 


A * D + E e ? '^ + E + E e"^ + E e - e - E, e ^ V 

P li o 1 2 3 ^ 5 6 

~ E e”^ - E e"^ ** 

7 8 


(22b) 


where E Q and D Q are the particular solutions. 


Proceeding along the same lines indicated in the proceeding paragraph, the 
integral equations for ^ and ^ ^ can be reduced to differential equations * 
Both of these differential equations have precisely the form of equation (21), 


with p ^ and ^ respectively replaced by ^ and |5 ^ or by p ^ and |6 2 ° • s °l uti °n8 
for this pair of simultaneous, linear, homogeneous differential equations are 


Pj =S *' 1 ' ( 3 3 =HePl 


( 23 ) 


Substituting these into the differential equations leads to a pair of linear, 

homogeneous algebraic equations which give eighth roots for p: + p , + p , + p , 

“ 1 - 2 “ 3 

+ p. and also that 
” h 

H t = G i , i = 1, 2, 3, U 

H ± *-Q ± , i * 5, 6, 7, 8 

With this information, the solutions for _ and . may be written as 

2 U 


-p,x 


-px 


-P X 


-P,x 


H n. 

^2 =G 3 8 + V + V 3 + V 




"V 

* G e A 4- Gi. e 


•PJ -P,Y -p Y 

2 -G ? e J = Gg e J 


(23a) 

(23b) 


P X pi 

The constants G , G , G^, and G are multipliers of terras such as e ^ , e , etc* 
1 2 h 6 

These exponentials approach infinity as X and Y approach infinity. Therefore, 
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to satisfy the condition that 


N and p 


are everywhere finite 


G 1= G 2=°5= G 6 =0 


( 21 .) 


To complete the solution, it only remains to determine the constants of 
integration which appear in equations (22a), (22b), (23a), and (23b). These may 
be found by returning to the four integral equations of the problem and intro- 
ducing the solutions for p p y and as given by equations (22a), (22b), 
(23a), and (23b). After the integrations are carried out, terms are grouped 
according to the exponential (e.g., e A , e c , etc.) which multiplies them. In 
this way, there are obtained lU linear algebraic equations for the lh unknown 
constants. These equations are solvable in closed form, and the final expressions 
are given in the Appendix. By inspection of these expressions, it is seen that 
the solutions depend upon the surface absorptivity o< and the inclination angle 
of the incident rays. 

With the solution for the radiosity B now at hand, consideration can be given 
to the energy absorbed in the cavity. The local rate of energy absorption per 
unit area is given by equation (U). Substituting the current definition of p from 
equation (16), the expression for q becomes 


.77 Ssi " J jf 


(25) 


The rate of energy absorption in the cavity as a whole can be obtained by integrating 
equation (25) along both walls of the cavity. Then, noting that the total rate 
of incoming energy is (S sin £ ) 1^, the apparent absorptivity can be evaluated 
from its definition (6) as 


c< 


Q. 


J&lQrj ¥ 
I ~o< 




fr 


Wc/Y 


( 26 ) 


where h/L^ has been replaced by tan £ according to equation (12). The solutions 
for as given by equations (22a), (22b), (23a), and (23b) may then be introduced. 
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and after the indicated integration is carried out, there is obtained 



Utilizing equation (26a) in conjunction with the expressions for p, D, E, 
and G as given in the Appendix, the apparent absorptivity has been evaluated for 
inclination an;les ^ ranging from 10* to 80° and for surface absorptivities 0< 
ranging from 0.3 to 0.9. These results are plotted in figure 5. From the figure, 
it is seen that the apparent absorptivity is greater as the inclination angle ^ 
decreases. This finding is made plausible by noting that the external radiation 
penetrates more deeply into the cavity as ^ decreases and must therefore undergo 
many diffuse reflections before appreciable amounts can escape. Thus, for ray 
hundles which .• rrive at deep diffuse cavities under small angles of inclination, 
there is. a vei/ appreciable cavity effect. As before, the cavity effect is 
accentuated vrth decreasing surface absorptivity. 

It is interesting to inquire as to whether the diffuse-surfaced cavity absorbs 
more effect! sly if the incoming energy is in a parallel ray bundle or is diffusely 
distributed. For deep cavities, this information can be obtained by comparing 
the curves <f figure 5 with the horizontal asymptotes of figure 2. From this 
comparison, it is found that rays which arrive at a small inclination angle ^ are 
more fullj absorbed than is diffuse incoming energy. The contrary is true for 
rays which arrive under a large inclination angle. Thewe findings may be illuminated 
by noting that a ray bundle arriving at a small inclination angle penetrates deeply 
into the cavilyj substantially deeper than the diffuse incoming stream. For a 
ray bundle arriving at a large inclination angle, the opposite is true. Further 
inspeauion of figures 2 and 5 reveals that when ^’.'£>30% the two typew of incoming 
radiation a:« absorbed about equally well. 
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The results of figure 5 were obtained from an analysis for the cavity of 
infinite depth. However, as already indicated, these results should also apply 
to finite cavities of sufficient depth. An approximate calculation was performed 
to determine the depth of finite cavity to which the infinite cavity results 
would apply to within 2%. These calculations were carried out by following groups 
of incoming rays as they diffusely reflected and re-reflected within the cavity. 
Account was kept of the energy absorbed at each surface contact and summation 
gave the total energy absorbed. For each inclination angle and surface absorp- 
tivity, the depth of cavity for which the apparent absorptivity is approximately 
98% of that given on figure 5 was found by trial and error. The cavity aspect 
ratios thus obtained are plotted in figure 7 as a function of inclination angle 
for parametric values of surface absorptivity. From the figure, it is seen that 
for cavities with surface absorptivities of 0.5 and larger, only moderate cavity 
depths are needed to achieve the absorbing power of a cavity of infinite depth. 

For cavities of lower absorptivity, deeper cavities are required, especially at 
small inclination angles. 

PARALLEL INCOMING RAIS, SPECULAR SURFACES 

The analysis for the case of parallel rays incident on a specular cavity 
may be carried through by a different and basically simpler approach than that 
given for diffuse surfaces. The derivation is facilitated by reference to figure 8. 
Here, there is shown a single ray entering the cavity, reflecting on the walls, 
and then emerging from the opening. A mirror image of the cavity as formed in 
the base surface is also shown at the left by the dashed lines. Suppose, for the 
moment, that the base surface is removed and that the ray can move freely through 
a cavity of length 2L. Then, the figure demonstrates that the reflection pattern 
in the dashed (image) portion of this cavity is identical to that which occurs 
in the actual cavity after the ray has been reflected from the base surface. 
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Thus, the number of surface contacts which a ray experiences within the actual 
cavity is equal to the number experienced in a cavity of length 2L open at both 
ends, plus one additional to account for reflection at the base. 

The ray which has been constructed in the figure has a special characteristic. 
Namely, that after making three contacts with the surface, it grazes the cavity 
wall as it leaves through the opening. Additional study of the figure indicates 
that rays which initially strike the surface in the region x<x will experience 
four surface contacts within the cavity, while rays which initially strike in the 
range x > x will experience three surface contacts. Thus, x divides the surface 
into two zones, such that rays striking within each zone have a uniform number 
of surface contacts. 

The number of surface contacts k experienced by rays which strike in the 

i 

region 0 x < x may be determined from the geometry of figure 8 as 

k = integral part £ 2 (L/L ± ) + 2 J (27) 

If the quantity in braces is precisely an integer, then k is equal to that integer 
minus one. For rays which strike the surface in the region x< x a L^, the number 
of surface contacts which are experienced within the cavity is k<= 1. Furthermore, 
from the geometry of the figure, it follows that the dividing point x is given by 

x/1^ = 2{h/\) - (k- 2) (28) 

When a ray experiences k contacts with the surface, the fraction of its 
energy which is thus absorbed is 

.j k-f v 

t- <*(l- o() -4- ol(i-ot) ■*- \ oUi-ct) (29a) 

This is a geometric series whose sum is 

1 •» (l~oO k (29b) 

With this, we can now calculate the apparent absorptivity c< from its definition, 

a 

equation (6). The energy absorbed in the cavity is 

S sinjif \jZ\l - (l~o() k ] + (L t =x) [l - (l-tf )*^! (30) 
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where S denotes the energy per unit area normal to the ray. The rate at which 
wnergy enters the cavity is S sin ^ 0 Then, introducing this information 
into equation (6) and taking x from equation (28), there is obtained 

o^ a = 2 - (k= 2)"[ j^l - (l~0<) k j + jk - 1 - 2 jl= (l-«) k " 1 j (31) 

It is easily verified that this expression is equally valid for L/l^ 1 as well 
as for L/L jl 1 1. 

Consideration may now be given to the computation of numerical results. 

Fran equation (27)* it is seen that for 0 £ L/l^— l/2, k = constant - 2j 

1/2 ^ L/L t 5'. 1, k ” constant = 3j 1< L/L^£=. 3/2, k - constant = U, and soforth. 

From this, it is clear that the number of surface contacts increases by one as a 

length increment 1/2 is added. Then, turning to liquation (31), it is seen 

that within a given range such as 0 ^ L/l ± ^ l/2, increases linearly with 

L/L o The same is true in the range 1/2 ^ I ,/L^ S 1, except that the slope of <X 

versus L/L is different (because k is different). Therefore, a graph of oi. 

x a 

is composed of a succession of straight line segment? which change slope at 
L/L ± = 1/2, 1, 3/2, 2, . . .o 

Figure 9 presents results for o( as a function of L/L. for parametric values 

a x 

of ck„ The previously-described changes in slope are clearly evident. From the 
abscissa, it is seen that L/h and ^ do not appear a» separate parameters, but 
rather, combine as the single parameter (h/L) tan J , Taking cognizance of this, 
it is evident from the figure that a given cavity (i.e., given L/h) absorbs energy 
more effectively when the rays arrive under a large angle of inclination (i.e., 
large )„ This is because such rays undergo a greater number of specular 
reflections than do rays which arrive at a small angle of inclination. This is 
precisely opposite to the characteristics of diffusely-” reflecting cavities, 
which are more effective absorbers for rays which arrive under a small angle of 
inclination (see figure 5). Additionally, for ray's arriving at a given angle of 
inclination (given / ), figure 9 indicates that the amount of energy absorbed 
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increases steadily as the cavity depth increases and ultimately approaches total 

absorption, i.e., o< * 1. For diffuse surfaces, it may be recalled that complete 

a 

absorption is not possible regardless of the depth of the cavity. 

It is interesting to inquire whether a specular cavity absorbs more effectively 
when the incoming radiation is a diffuse stream or a parallel ray bundle. From 
a study of figures 2 and 9, it is found that parallel rays arriving at a large 
angle of inclination are more strongly absorbed than is diffuse incoming energy. 

On the other hand, rays arriving under a small angle of inclination are less 
strongly absorbed than is diffuse energy. 
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APPENDIX 


Expressions for p, D, E, and G. 



Then D = Z 2 E = Z /«o_ 
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where i, j, and m are given by the following table: 

k i i m 

112 1 
2 2 11 

5 3 h 2 

6 1 : 3 2 
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jf, and m are given by the following table* 
k i i m 

3 12 1 

1 * 2 1 1 

7 3 U 2 

8 U 3 2 
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Fig. 5 Apparent Absorptivity for Incoming Parallel Rays in a Diffuse Cavity 
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Fig. 7 Aspect Ratio of Finite Cavity at Which Infinite Cavity 
Results Apply Within 2 Per Cent 




1 


I 


Fig. 8 Incoming Parallel Rays in 
Cav ity 




a Specular 



Fig. 2 Apparent Absorptivity for Diffuse Energy Entering a Diffuse Cavity 




1,2,3: SPECULAR 4: Dl FFUSE ( BLACK ) 


Fig. 3 Images Formed by Diffuse Incoming Radiation in 

a Specular Cavity 




2 3 4 




<X=D.I 



Fig. 4 Apparent Absorptivity for Diffuse Energy Entering a Specular Cavity 





Fig. 5 Apparent Absorptivity for Incoming Parallel Rays in a Diffuse Cavity 






y ( Degrees ) 

Fig. 7 Aspect Ratio of Finite Cavity at Which Infinite Cavity 
Results Apply Within 2 Per Cent 




Fig. 8 Incoming Parallel Rays in a Specular 
Cav ity 




